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Abstract 

We compare the resonant and non-resonant contributions in various regions of phase 
space for the D — > Kit semileptonic transition amplitude, computed in a chiral model 
which incorporates the heavy quark symmetry. Remarks on the significance for experi- 
ment and for chiral perturbation theory are made. 
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1. INTRODUCTION 



In this note we shall examine the weak hadronic current matrix element for the decay 
D° — > K~7T°e + i> e using a chiral Lagrangian which incorporates the heavy quark symmetry 
Previous works have treated this process in a chiral model which includes only light 
pseudoscalars [1] and in chiral models with both light pseudoscalars and light vectors 
present but with the approximation that the decay be replaced by D° — > K*~e + v e [2-6]. 
Here we will consider both contributions together. This is interesting, of course, in its own 
right. It also holds some interest for the question of what is the best way to incorporate 
vector mesons in the chiral perturbation theory program [7]. It may actually be easier to 
investigate this question in the framework of chiral light-heavy interactions rather than 
light-light interactions since the "heavy end" might eventually be under better control. 
For the present process we find that there is no region of phase space in which the light 
vector K* piece does not make a non-negligible contribution and that it is typically very 
dominant. This perhaps suggests the adoption of a framework in which light pseudoscalars 
and light vectors are treated together from the beginning. 

The detailed points include the discussion of the way the chiral theorem is maintained 
in the appropriate unphysical limit and the treatment of the phase space kinematics for 
the hadronic matrix element. The transition amplitude is given in Section 2 (see also 
the Appendix) after our notation has been introduced. Section 3 contains the kinematics 
and the comparison of the resonant and non-resonant pieces of the amplitude in various 
regions of phase space. Some caveats and remarks on the experimental aspect of K* 
dominance are given in Section 4. 
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2. Decay Amplitudes 



We are interested in the chiral invariant interactions of both the light pseudoscalar 
nonet <p(x) and the light vector nonet p M (x) with the heavy meson field. We shall follow 
the notation of Ref. [3]; other treatments include Refs. [2, 4, 5, 6]. The chiral interactions 
involving only the light pseudoscalars were discussed in Ref. [8]. Using the "heavy field" 
H(x) which contains both heavy pseudoscalar as well as heavy vector pieces, the leading 
order (in heavy meson mass M and in number of derivatives of the light fields) strong 
interaction is compactly written as [3]: 

-^Aight - heavy = iV„ Tr {H[d^ - iagp^ - i(l - a)v^]H} 

+ id Tr [H lfll5PlI H] + — Tr [H^^F^H], (2.1) 

ui v 

wherein m v is the light vector meson mass introduced just to keep the coupling constant 
c dimensionless and 

WM^^i^O, (2-2) 

with the chiral matrix £ = exp(i<p/ F n ). Furthermore our normalization convention sets 
F w ~ 132 MeV. is the heavy meson 4- velocity and the heavy field H here is taken 
to have the canonical dimension of one. g ~ 3.93 is the light vector-light pseudoscalars 
coupling constant. Heavy quark symmetry breaking terms, SU(3) symmetry breaking 
terms as well as the chiral Lagrangian of the light sector have not been explicitly written. 

Notice that the light-heavy interaction (2.1) is characterized by the three dimensionless 
coupling constants a, c and d (denoted g in [8]). The choice a = 1 corresponds to a 
natural notion of light vector meson dominance. This choice sets to zero the coefficient of 
v^ — i/ (2F^)(0(9 M — <9 M 00) + ... so that two pseudoscalars in a p-wave state can only be 
emitted through an intermediate light vector particle from a single heavy meson vertex. 
Whether, in fact, a« 1 remains to be determined. 
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For our present application we also require the four fermion effective weak interaction: 

'-w — A* n ' 
J ( ~ ] =iv e l^(l +l 5 )e + ... , 

4 +) =tV; s s^(l + l5 )c+..., (2.3) 

with usual conventions [3] and where V cs is the Kobayashi-Maskawa matrix element. The 
chiral covariant realization of the left handed hadronic current, in terms of heavy 
and light meson fields is 

4 +) / V cs = F D [d^D b + ia'~gD aPimb + i(l - a')D a v mb + MD; b ](£% + ... , (2.4) 

where the SU(3) triplet fields (D 1 , D 2 , D 3 ) stand for D°,D + Df) and similarly for the 
heavy vectors. (Eq. (2.4) is the same as (4.6) of [3], but we have redefined (a + a') by 
a' to avoid confusion). In (2.4) a' is a new dimensionless coupling constant (which scales 
however as M) characterizing the phenomenological hadron weak current. We can rewrite 
(2.4) in the heavy quark limit as 

tfPlK = Tr [ 7 „(1 + l5 )H a ](^U + \FdcI Tr ( l5 H a )(~g Pfmb - v, ab )(^) b3 + ... . 

(2.5) 

Now let us compute the hadronic matrix element for the process D°(p) — > K~(p') + 
n°{p") + e + (g e ) + v e {q v ) We define the 4-momentum 

q = q e + qu=p-p' - p" , (2.6) 
and employ the following form factor decomposition: 

^ Po p' p^(p')Apn\Jl + yv;jD°(p)} 

= ~i[q»r + (p' + p")^ + + (p' - p")^ + he^Ppp'M . (2.7) 
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The first term on the right hand side of (2.7) does not contribute to the net weak amplitude 
since dotted into the leptonic factor vanishes for zero lepton masses. Hence the form 
factor r is irrelevant. There are contributions to the form factors both with and without 
intermediate light vector meson K* poles. The "non-resonant" (NR) diagrams without 
the intermediate K* are discussed in the Appendix. Taking the leading order in M 
contribution to each form factor yields 

/ F D \ dM 



K) 



NR — 



2y/2F%) A-V-p»' 



{uj - )nr = (^)(-A^7 + a/ )' 

F od 2 1 1 

NR ~ V2F* A*-V- (P + p") A-V-p»> (2 ' 8) 

where V M = p^/M is the 4-momentum of the initial D°, A = M(D*) - M(D) and 

A* = M(D* S ) - M(D). uo + and both scale as M 1/2 while h (since it multiplies p p ) 

scales as M" 1 / 2 . Eq. (2.8) differs from an earlier calculation [1] in the model with light 

pseudoscalars only by the a' term in (u_)nr- ol' is due to the presence of light vectors, 

as may be seen from (2.4). Chiral covariance demands an additional pseudoscalar piece 

when we add the light vectors. The analogous "strong" parameter a does not appear in 

(2.8) but does show up in the form factor r NR given in (A7). 

The computation of the diagrams containing K* poles can be simplified by making 

use of earlier results on the D — > K* weak current matrix element [2-6]. From section 5 

of [3] we obtain the leading large M contribution in this case as 



^p~^ (K*-(k,e)\J^/V; s \D°(p)} 



J^ko(K* (k,e)\JZ/V*\D + (p)) 



iF D e u 



2cM Vpk c agV u q p 

« 9<J^u H 



A*-V-k A s -V-k 



(2.9) 



wherein % = p^ — k^. Additional corrections due to higher derivative interactions [3,4], 
loops [4] and excited heavy states [2,6] have been discussed in the literature but (2.9) 
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seems sufficient for our present purpose. The pieces of the D°(p) — > K~(p') + ir°(p") 
transition matrix element in (2.7) involving K* poles can be found from (2.9) simply by 
defining 

K=p', + p'1 (2-10) 

and replacing e u by the factor, 

g ( m v \ 2 [p' u - - fc^(m^ - ml)/m 2 K *\ 
v 2V2\F n g) k 2 + m 2 K *-im* K r K * ' { ^ ii> 

wherein the combination {jr~) 2 is numerically close to 2.0. Eq. (2.11) is the product of 

the vector -2 pseudoscalar coupling constant and the K* propagator. Since the K* can go 

"on shell" we include the conventional width correction in the denominator to maintain 

unitarity. Decomposing the "resonant" K* pole amplitudes (subscript R) into the form 

factors defined in (2.7) gives: 

<<*+>« = - ± ^ L <"->«. 

( \ f F D \ a! ml 



\2^/2F%) k 2 + m 2 K *-im K *T K * ) 

= 2c~gF D /mA 2 1 1 

R V2m v {F^gJ A* - V ■ (p> + p") k 2 + m 2 K . - im K »T K , 1 ' ' 

The net results for the D° — > K~n° current transition form factors in the limit of 

heavy charmed particles interacting with "soft" light pseudoscalars and vectors are given 

by the sums of the corresponding terms appearing in (2.8) and in (2.12). The precise range 

of validity of the concepts of a soft light pseudoscalar and (especially) a soft light vector 

are probably best left to a comparison with experiment. Certain results, which depend 

only on the existence of spontaneously broken chiral symmetry, are designated "current 

algebra theorems" and should hold for zero mass pseudoscalars as their 4-momenta go to 

zero. In particular we may note that (uj+)r vanishes in this limit and that {oj~)r (when 

we set m 2 = m 2 K » and IV* = 0, the latter corresponding to a pure effective Lagrangian 

computation which thereby satisfies chiral symmetry by construction) cancels off the a' 
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piece of (u-)nr- In this unphysical limit u + = —uj- [9]. We should stress that the fact 
that a' cancels out in this limit does not mean that it is not an important parameter for 
describing the decay; actually, it turns out to be the most important parameter. 

To proceed with a comparison of the resonant and non-resonant contributions we need 
at least a rough idea of the magnitudes of the "strong" parameters a, c and d as well as the 
"weak" parameter a'. Bounds on the value of d have been obtained in the literature [10] 
which agree with a simple estimate [3] based on pole dominance of the D — > K transition 
form factor: 

d^0.53. (2.13) 

Information about c and a' may be obtained by comparing (2.9) with experimental infor- 
mation on the decays D — > K*e + u e . Eq. (2.9) is expected to be most reliable for "soft" 
fT*'s, which implies that (— q 2 ) should be as large as kinematically possible. Now the ex- 
perimental data is analyzed in terms of form factors characterizing (2.9) which have the q 2 
dependence, M* 2 /(M* 2 + q 2 ). Evaluating (2.9) at the extreme value -q 2 — (M — m K *) 2 , 
assuming the above q 2 dependence to extrapolate to small — q 2 and comparing with the 
experimental values [11] yields the estimates 

\a'\ « 1.73, |c|«1.6. (2.14) 

a is not determined from this analysis. However a study of the binding energy of heavy 
baryons as solitons of the meson Lagrangian involving (2.1) shows that (5.3) of Ref. [12] 
may be fit with the choice 

= 0.53, c=1.6, a = -2. (2.15) 

We stress that the numerical estimates (2.13)-(2.15) are preliminary in nature; however 
they should suffice to draw qualitative conclusions. 
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3. Comparison of Light Pseudoscalar and Light Vector Amplitudes 



First we must discuss the kinematics associated with the hadronic matrix element 
(2.7). For the decay D°(p) — > K~{p') + n°(p") + e + (q e ) +v e (q v ) there are five independent 
dynamical variables needed to specify the momentum configuration; a recent discussion 
is given in Ref. [1] (see also [13]) based on earlier treatments [14] of k e4 decays. It is 
somewhat simpler if we confine our attention just to the hadronic part (2.7). Then, for 
each value of the invariant lepton squared mass — q 2 we have in effect a three body final 
state. The classic Dalitz plot analysis shows that the kaon energy E' and the pion energy 
E" (both in the D° rest frame) are sufficient. Altogether this gives E', E" and q 2 as a 
possible complete set of variables to describe (2.7). Since k 2 plays an important role in 
the K* pole amplitude we shall choose the alternative set 

(k 2 ,q 2 ,E"). (3.1) 

A convenient formula relating the two sets is 

M 2 - k 2 4- a 2 

E' + E"= M 2 k M + q =Q(k 2 ,q 2 ), (3.2) 
where M is the D° mass. With the neglect of the lepton masses, —k 2 satisfies 

(m n + m k ) 2 < —k 2 < M 2 . (3.3) 

For a given k 2 , q 2 must satisfy 

< —q 2 < (M — v / — A^") 2 • (3.4) 

The region defined by (3.3) and (3.4) is illustrated in Fig. 1. Now, for each value of q 2 
there is a Dalitz plot boundary in the k 2 — E" plane obtained by rewriting the condition 
|p' ' P"l — Ip'IIp"! i n terms of the set (3.1). The allowed regions, corresponding to 
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horizontal slices perpendicular to the page in Fig. 1, are illustrated in Fig. 2. A simple 
formula can be given for the "average" value of E" at each (k 2 , q 2 ): 



E"(k 2 ,q 2 ) = \\El^k 2 , q 2 ) + E'Uk 2 ,q 2 ) 



2 



1 + 



mir — m 



K ~ >""K 

k 2 



(3.5) 



Notice that Fig. 1 does not correspond to a fixed E" slice but rather represents the 
projection of the (k 2 ,q 2 ,E") boundary surface into the (k 2 ,q 2 ) plane. 

The expressions for the form factors in the heavy quark limit given by (2.8) + (2.12) 
are expected to be most reliable for large (— q 2 ), corresponding to "soft" n° and K~ 
particles as well as, in the spirit of the present Lagrangian, soft K* particles. Refering 
to Fig. 1, we see that the soft K* condition suggests that we consider the effective 
Lagrangian expressions to hold in the region where — q 2 is greater than around 0.5 GeV 2 . 
This represents a fairly large portion of the entire phase space. Nevertheless it is not 
unreasonable from an experimental point of view inasmuch as the Particle Data Group 
tables state [15] that "it is generally agreed that the Kire + v e decays of the and D° are 
dominantly K*e + v e . n Experimentally, it is also known that the amplitudes are damped 
for decreasing (— q 2 ); this can be seen from Fig. 1 to decrease the importance of the larger 
— k 2 region. 

From a theoretical point of view it is very interesting to compare the pseudoscalar 
and vector contributions to the hadronic form factors. Nature tells us, of course, that the 
vector contributions dominate when one folds in the leptons and integrates their "squares" 
over all phase space. Nevertheless, it is important to get an idea of the "local" ratios of 
the two contributions in various kinematic regions. The phenomenological analysis of 
Ref. [16], for example, shows that it is the form factor proportional to e M in (2.9) (i.e, 
the a' term) which mainly supplies the total width for D° — > K*~e + V e . This form factor 
contributes to both (cu+)r and (uj^)r in (2.12) but (w-)r is clearly the dominat one; 



(uj+)r would vanish when the k and the 7r masses are equal or neglected. Hence we will 
focus our attention on the oj- form factor. It is convenient to rewrite it as follows (in the 
D° rest frame): 



0J_ 



'V 



a' 1 



m 



K* 



dM 



(3.6) 



2V2F%[ V /c 2 + m 2 K * - im K *T K * ) A + E" _ 

We noted earlier that the coefficient of the a' term (which corresponds to the extra piece 
introduced by adding light vectors to the effective Lagrangian) vanishes in the unphysical 
— > limit (with IV* = 0) in agreement with expectations. Let us then consider the 
ratio of the magnitude of the a' term to the magnitude of the last term: 



ratio 



;a + e") 

M 



(^ 



m 



K* , 



m 2 T 



1/2 



(3.7) 



Only a' and d are not very well known; we will use the estimates in (2.13) and (2.14) 
which are certainly qualitatively reasonable. There is a lower bound for the above, 



ratio > 



A + 



0.5, 



{3.1 



M m? K * — [m K + m n ) 2 

where Y K * was set to zero for simplicity since it has a negligible effect when (—A; 2 ) is 
as small as possible. This result indicates that there is no region in which the vector 
contribution is negligible compared to the pseudoscalar contribution. So if one were to 
make the usual derivative expansion of the chiral perturbation theory approach, there 
would be large corrections to the first order results due to the existence of the K*. Of 
course, the pseudoscalar piece takes on its largest value near the cusp in Fig. 1. 

Since (3.7) depends only on k 2 and E" it is convenient to display curves of constant 
ratio in Fig. 2. We notice that there is only a very small region for which the ratio is less 
than unity. Because of the small width of the K* (50 MeV) the ratio rises dramatically 
to around 15 at — k 2 = m 2 K *. It is seen that the large ratio of amplitudes persists over a 
non-negligible region. 
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The ratio 



of the magnitudes of the vector and pseudoscalar contributions to the 



h-NR 

weak vector current form factor is seen from (2.8) and (2.12) to be the same as the ratio 
(3.7) when we multiply the latter by the factor 

mK * M (3.9) 



a'g 



(A; 4 + m 2 K S 2 K ,) l l 2 ' 

This is numerically around 1.0 in the phase space region of interest so the K* contribution 
is also dominant for this form factor. 

Tu sum up, in the large (— q 2 ) region (optimistically as large as — q 2 > 0.5 GeV 2 ) the 
K* contribution overwhelms the pseudoscalar contribution nearly everywhere. Even for 
the very largest (— q 2 ), where the soft pseudoscalar results are expected to be most signif- 
icant, the K* amplitudes are relatively sizeable. This situation would seem to suggest the 
desirability of a modified chiral perturbation theory program in which both pseudoscalars 
and vectors are retained in the effective Lagrangian from the very beginning. Some re- 
cent discussion of this point of view has been given in Ref. [17]. In the case where one 
is considering a non-strange transition matrix element (like B — > mi) we should replace 
in (3.8), {rriK + rn^) by 277^ and m^.by m p . Then there would be a very small region 
in which the pseudoscalar piece might be considered dominant, but the overall picture 
would be qualitatively similar. 

4. Remarks 

1. Of course, it would be a better approximation to deal with the B rather than the 

D meson as an example of a heavy field. Similarly it would be a better approximation to 
restrict the light particles to the non-strange ones. Thus a similar analysis (with basically 
identical formulas) would be somewhat cleaner for B — > ixix current transitions rather 
than D — > Kn transitions. Nevertheless there is at present much more experimental data 
for the D — > Kit case and it is quite identeresting in its own right. 
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2. As noted after (2.9), the formula we are using for the resonant contribution ne- 
glects a number of effects which are subleading from the light meson point of view but 
may be necessary to take into account if the relatively small form factor which would be 
proportional to in (2.9) is actually non-zero. However we expect that our approxi- 
mation is sufficient for the purpose of comparing the relative strengths of the vector and 
pseudoscalar contributions. 

3. The ratio of the non-resonant part of T(D — > Kne + i> e ) to Y(D — > K*e + u e ) is of 
evident experimental interest. The precise meaning of this quantity depends on the man- 
ner in which it is extracted from experiment. Apparently, there is no universal method. 
The most straightforward way is simply to define the resonant contribution as everything 
within a certain band of — k 2 surrounding m 2 K *. This definition has, however, the mis- 
leading feature that is counts non- resonant background near the peak as resonant. It is 
particularly easy to apply this definition to the present case when one makes the reason- 
able approximation that the entire amplitude is dominated by the K* pole diagrams. As 
recently illustrated for a different decay in Ref. [13] both the phase space and the squared 
amplitudes factorize in this approximation so that one obtains 

r(D° -> K-*°e+v e ) « r( ^" T T K ~^ ] T(D° -> K*~e^ e )- f°° (4.1) 

1 K* TT J-oo X A + 1 

where x is defined by k 2 + m 2 K , = xY^ttik*- If the "resonant" region is taken to be 
that range of —k 2 for which \\/—k 2 — m K *\ < NT K * then integrating (4.1) yields the 



non-resonant /resonant ratio to be about 



TT 



1. (4.2) 



2 tan-\2N) 

This gives a 19% non resonant contribution for = 2 and 12% for = 3. The extent 
to which the data obeys (4.2) as a function of might be considered a measure of how 
good is the K* dominance. The fact that we do get a non-resonant contribution at all 
with pole dominance is, of course, an artifact of its definition. However, it also illustrates 
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the difficulty in giving a meaningful experimental definition of the non- resonant /resonant 
ratio. Since the "theoretical" resonant and non-resonant amplitudes have the same order 
of magnitude outside the resonance region, the above numbers may give a satisfactory 
rough order of magnitude estimate for any reasonable definition of the ratio. 
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Appendix 



For the process D° — > K 7r°e + z/ e we require the hadronic matrix element 



8p M(K-(p')7r°(p")\J^V: s \D (p)), (Al) 

where the states are normalized in a unit volume. We shall proceed by first using the 
ordinary, rather than "heavy", meson fields to compute (Al) (see section V of Ref.[3]) and 
then take the heavy quark limit. The contribution to (Al) from the contact (non-pole) 
diagram is: 

^^\p, + (a'-l)(p^-p'p}. (A2) 

The D+ pole diagram, wherein D°(p) -> K~{pf) + n°{p") + D+(q) followed by D+(q) -> 
e + z/ e , contributes the term 

iF D (l-a) (p - p") •jp + q) ( aq\ 

2V2F* M 2 + q 2 % ' KA6) 

in which M s is the Df mass and q^ = p^ — (p' + The D*° pole diagram, wherein 

D° — > 7T° + _D*° at the strong vertex followed by -D* — > X - + at the weak vertex, 

contributes: 

-iV2M 2 F D d [pl + {p-p")^p-p")-p"/M* 2 ] 

F 2 ( p _ p //)2 + M *2 ' y ) 

in which M* is the -D* mass. There are also two double-pole diagrams. The first features 

D° — > 7T° + -D* at a strong vertex followed by D*° — > fT~ + _D* + at another strong vertex 

which, in turn, is followed by D* + — > e + z/ e at the weak vertex. The contribution to (Al) 

is 

-i2V2d 2 M 2 F D ^ ap vP P p' u p'L (A) 

F 2 [Mf + q 2 }[M* 2 + (p - p") 2 Y 1 ' 

in which M* is the D* + mass. Finally, the diagram with D° — > 7r° + -D* , -D* — »■ + 
Z}+, — > e + z/ e gives 

-i2V2F D M 2 d 2 q^p' ■ p" +p'-(p- p")p" ■ (p - p")/M* 2 } 



F 2 [q 2 + M 2 } [{p - p") 2 + M* 2 ] 
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(A6) 



In order to obtain a model independent result corresponding to M — > oo we should delete 
terms which fall off faster than M 1//2 , the scale behavior of (Al). For example, using 
F D ~ M- 1 / 2 , and a' ~ M, (A2) -> 



-iF D 



[p^ + a'ip'-p"),]. {A3!) 



In taking the limit of (A3) we set p M = MV^ and throw away terms of quadratic order 
in p' and p". The resonance denominator of (A3) becomes 2M[A S — V • (p' + p")\ with 
A s = M s — M and (A3) -> 

2V2F 2 A s — V • (p' + p") { } 

The other diagrams are treated similarly. With the form factor decomposition defined 
in (2.7) we obtain the results listed in (2.8) for the three experimentally significant form 
factors. For the sake of completeness we also give here: 

F D 



r N R = 



1 (1 - a) V-(p'- p") dV ■ p" d 2 (p' ■ p" + V ■ p'V ■ p") 



2 2F2 A s -V-(p' + p") A - V ■ p" [A s - V ■ (p' + p")](A - V ■ p") 

(A7) 
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Figure Captions 

Fig. 1 Projection of the three dimensional (k 2 ,q 2 ,E") phase space boundary into the 
(k 2 , q 2 ) plane. 

Fig. 2 Phase space boundaries in the k 2 — E" plane at various values of q 2 . In increasing 
order of size the closed curves correspond respectively to — q 2 = 0.925, 0.725, 0.525, 
0.325 and 0.025 GeV 2 . Also shown are points on the contour lines on which the ratio 
in (3.7) takes on fixed values; the circles, crosses and squares correspond respectively 
to the ratio equal to 1,3 and 9. 
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